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THIN  CYLINDRICAL  SHELLS  UNDER  LOCAL  AXIAL  LOADINGS 

ABSTRACT 


The  linear  equations  due  to  Fl'ugge  for  thin  elastic  circular  cylin¬ 
drical  shells  are  solved  by  use  of  the  finite  Fourier  transform  for  the 
case  of  distributed  tangential  (axial)  loading  with  simply  supported  edge 
conditions.  This  solution  supplements  prior  results  obtained  by  Bijlaard 
for  radial  (normal)  distributed  loadings. 

Combining  Bijlaard1 s  results  with  the  solutions  obtained  in  this  report 
permits  discussion  of  localized  loadings  having  components  both  normal  and 
tangential  to  the  shell’s  middle  surface.  Such  an  application  is  discussed 
and  numerical  results  are  provided. 
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INTRODUCTION 


It  frequently  happens  in  pressure  vessel  and  missile  construction 
that  an  axial  load  is  transmitted  into  the  wall  of  a  thin  cylindrical 
shell  by  means  of  a  heavy  ring  attached  to  the  cylinder.  This  ring  dis¬ 
tributes  the  load  fairly  uniformly  into  the  cylinder.  For  reasons  of 
economy  or  weight-saving,  it  is  sometimes  desirable  to  replace  the  solid 
ring  by  a  series  of  brackets  or  angle  irons.  The  resulting  nonuniform 
distribution  of  axial  load  can  result  in  a  significant  increase  in  the 
local  stress  level  in  the  cylinder.  This  situation  thus  provides  an  inter¬ 
esting  and  useful  problem  in  the  asymmetric  deformation  of  thin  shells. 
Solutions  and  numerical  results  for  similar  localized  radial  (normal)  and 
tangential  (in  the  circumferential  direction)  loadings  on  simply  supported 
cylinders  have  been  presented  by  Bijlaard  Cl, 2].  However,  the  related 
problem  of  a  distributed  tangential  loading  in  the  axial  direction  was  not 
discussed.  It  is  our  intention,  then,  to  supplement  these  results  by 
examining  the  solution  and  presenting  some  numerical  results  for  localized 
axial  loadings  on  thin  cylindrical  shells. 

In  the  specific  application  described  earlier,  the  axial  load,  is  not 
applied  directly  to  the  shell's  middle  surface  but  usually  at  a  distance 
greater  than  half  the  thickness  away  from  the  middle  surface.  Thus,  there 
is  also  developed  a  radial  load  with  an  axial  variation,  i.e.,  a  longi¬ 
tudinal  moment,  to  use  Bijlaard's  nomenclature.  The  desired  solution  to 
the  physical  problem  is  then  a  properly  weighted  combination  of  these  two 
loading  conditions. 

In  the  ensuing  discussion,  the  governing  differential  equations  of 
equilibrium  selected  are  those  due  to  Flugge  rather  than  those  employed  by 
Bijlaard.  This  approach  affords,  as  a  by-product,  an  opportunity  to  com¬ 
pare  results  of  Flugge 's  nonhomogeneous  equations  with  those  of  Bijlaard's 
for  a  wide  range  of  shell  parameters.  Hoff  [3]  and  Kempner  [4]  have  made 
similar  comparisons  with  Donnell's  equations  for  the  corresponding  homo¬ 
geneous  equations  and  for  line  loadings. 
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Recently,  Nash  and  Bridgland  [5]  discussed  line  radial  loadings  on 
thin  cylindrical  shells  using  Flugge's  equations  and  suggested  the  use 
of  the  finite  Fourier  transform  in  obtaining  solutions.  This  technique 
has  been  employed  in  the  present  discussion  of  distributed  axial  loads. 

BASIC  EQUATIONS  AND  FORMULATION  OF  PROBLEM 

We  consider  a  thin  circular  cylindrical  shell  of  radius  a  and 
thickness  h.  Young's  modulus  and  Poisson's  ratio  are  denoted  by  E  and  v, 
respectively,  while  x*,  0,  and  Z  denote  coordinates  in  th<?  axial,  circum¬ 
ferential,  and  radial  (positive  when  inward)  directions.  Corresponding 
displacement  components  for  a  point  on  the  shell's  middle  surface  are 
u*,  v*,  and  w*.  It  is  convenient  to  deal  with  the  dimensionless  dis¬ 
placement  components  u,  v,  and  w,  defined  by  u  =u*/a,  vevx/a,  and  w=w*/a, 
as  well  as  with  dimensionless  middle -surface  coordinates  x=x*/u  and  9. 

Fldgge's  equations  for  the  equilibrium  of  thin  circular  cylindrical 
shells  [6]  may  then  be  written  in  the  dimensionless  form 

u"  +  (^)V  +  (^)  v'-uw'  +  k2  [w"  +  (^)  (u  -  »')] 

V  +  (^2)  v  +  (^)  u'  -  w  +  k2  [(^2)  w”  +  V"] 

vH;  +  2w  +  w  +  (•2~)  v  +  w'"  -  (-i^)  u'  +  —  [w-v-uu'] 

2  2  k3 

where  dots  indicate  differentiation  with  respect  to  0,  and  primes  indicate 
differentiation  with  respect  to  x;  p  ,  pQ,and  p  are  the  components  of 
applied  surface  loading  in  the  axial,  circumferential, and  radial  direc¬ 
tions,  k2  =  h2/l2a2,  C  =  Eh/(l-u2),  D  =  Eh3/l2(l-y3),  and 
vSj  *  y ""  +  2  w"  +  w*  . 

When  Equations  1  have  been  solved  for  u,  v,  and  w  for  a  prescribed 
loading  condition,  stress  resultants  and  moments  in  the  cylinder  may  be 
obtained  from  the  following  relations: 


aP, 


a2p. 


(1) 


asvx 

~D~ 
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N  =  C  [u'  +  u  (Vw)]  +  D  w" 

DC 

Nq  =  c  [v+  wu'-w]  -  D  [w  +  w] 

Nxe  =  C  [u  +  v*]  +  D  (^)[y*  +  w'] 

Nex  =  C  (~^)  [u  +  v]  +  D  (^p)[u  -  w  ] 

(2) 

M  =  -Da  [w"  +  v  +  u1  +  v  v] 

DC 

M  =  -Da  [w  +  w  +  v  w"] 

0 

Mx0  =  "®a  (1"UK*'  +  v1] 

MQx  =  -Da  (l-„)0*  +  |  (V-u)] 

D  =  D/a2. 

We  consider  the  following  specific  loading  and  boundary  conditions. 

An  axial  load  of  intensity  p  is  distributed  over  a  set  of  N  rectangular 

DC 

areas  (pads)  periodically  located  around  the  circumference  of  the  cylinder 
(see  Figure  l).  Let  the  length  in  the  axial  direction  of  each  pad  be  2a6 
and  the  x  coordinate  of  the  center  of  each  pad  be  x^;  let  the  length  of 
each  pad  in  the  8  direction  be  2a 8i ,  while  the  0  coordinate  of  the  center 
of  the  first  pad  is  0=0.  The  distributed  applied  load  is  assumed  to  be 
balanced  by  a  uniform  axial  compression  N  at  the  base  of  the  cylinder. 

DC 

This  desired  loading  system  may  be  considered  as  the  superposition  of  the 
following  two  systems: 

1.  a  rotationally  symmetric  distribution  of  uniform  axial  load 
acting  over  a  band  of  length  2a6,  centered  about  x=Xq,  and  balanced  by  a 
uniform  compression  N  at  the  base  of  the  cylinder. 

DC 

2.  a  self-equilibrating  distribution  of  tangential  axial  load  acting 
over  the  same  area  (see  Figure  2). 

The  solution  to  the  first  of  these  systems  is  trivial.  Hence,  the 
given  problem  is  reduced  to  the  solution  of  equilibrium  equations  when 
Pq=P^=0  and  px(x,8)  has  the  form  of  a  step  function  in  x  for  0=0q  and  a 
rectangular  wave  function  in  0  for  xq-6  5  x  £  xq  +  6  (see  Figure  3)* 


\ 
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Algebraically,  then, 

Px  (M)  =  Pi  2(mJ.)(ex+qa)  -  6i  2  ©  s  2(m-l)(0l+Qa)  +  0i 

(3) 

=  -Ps  2(^1)  (01+08 )  +  01  ^  ^  ^(Si+eg)  -  0i 

where  m=l,2,...  N,  0i+02  =  rt/N,  xq-6  sx  ^  xq  +  6. 

For  x  outside  this  range,  p  vanishes  identically.  The  condition  that 

A 

the  load  he  self -equilibrating  is  that  px  0X  =  ps 02 .  Finally,  we  consider 
boundary  conditions  of  simple  support  at  the  ends,  i.e.,  w  =  v  =  W  =M  =0 
at  x  =  0,  L/a. 

METHOD  OF  SOLUTION 


We  assume  the  following  forms  for  the  displacements: 


vs  2  Fjn  (9)  sin  (kx.) 
n=l 

v  =  2  Fan  (0)  sin  (Xx) 
n=l 


u  =  F3o  (0)  +  2  F3n  (0)  cos  (Xx) 
n=l 


where  X  = 


rota 

L 


This  form  ensures  simply  supported  edges  at  x=0,  L/a.  We  further  expand 
the  applied  load  pv  (x,0)  in  a  Fourier  cosine  series  in  x,  obtaining: 

aj0) 

+  2  t. 

n 


x 


Px  (*>e)  = 


+  2  a„  (®)  cos  (kx) 
n=l 


an  (0)  =  g-  P  (0)  cos  (Xxo)  sin  (X§)  (5) 

a  (e)  =  itfiSELSl  . 
o  '  '  L 
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Insertion  of  these  expansions  into  the  equilibrium  equations  (l)  leads  to 
a  set  of  3n  ordinary  differential  equations  for  F^(  6)  whose  solution  may 
rather  compactly  be  obtained  by  means  of  the  finite  Fourier  transform. 

The  essential  facts  regarding  this  transform  may  be  listed  as  follows. 

The  finite  Fourier  transform  of  a  function  f ( ©)  is  defined  as: 

f*(s)  =  F  [f(e)3  =  ~  f  f(e)  e"iSe  de. 

0 

The  inverse  transform  is  given  by 

f(e)  =  F”1  [f*(S)3  =  E  f*(S)  eiS9  , 

S=-“ 

while  the  following  useful  relation  holds  [7] 

F  [§  ]  =  iS  f* (S) . 


Application  of  the  transform  then  reduces  the  equilibrium  equations 

to  a  set  of  3n  algebraic  equations  to  be  solved  for  the  transforms  of 

F(e): 
mn 


Pm  (s) 


k3q4  +  2k3q3(l-\3)  +  [l+ka(l+A4)]j  +  Fsn*(S)  ■ -q[l+k8\s(^)]j 


+  F3/(S)  Ak3q3(^)  +  1  (v+*3k3) 


=  0 


Pm  (S) 


-q  [l+A3k3 (rTr) ]|  +  P8n*(S) 


+  j- 


9»-ls  (^)(l+3  *s) 


*  <T*» 


=  0 


Pm  (S) 


I2**2  (AH)  +  i  (, hW)}  +  Fan’(S)  Uq  (SS) 


2 


+  F3n*(S)|x3- 


q3  (l+k3)(^)l  = 


Pp(e) 


(6) 


where  q=iS. 


The  transform  of  p(8),  the  8-dependence  of  the  load  function  is  easily 
found  to  be 


F  [p(8)]  =4  L  sin  (S  0,)  Jx  *s 


•  fa  a  \  +Q3 )  isi 

-p2  sm  (S  0S)  e  '  / \  i  Y3 /  j 


(7) 


f 


=  0 


s  =  kN 


k=l,2,3* • • 
s  £  kN. 


With  this  information,  simple  algebra  permits  an  explicit  solution  for  the 
transforms  of  F  (8).  Inserting  these  into  the  inversion  formula  and 
taking  account  of  the  odd  or  even  character  of  the  resulting  expressions 
leads  directly  to  the  following  results  for  3^(0) ; 


F  (0)  = 
mn'  ' 


Fan(6)  = 


2a(anjW(p1+ps) 

Crtki 


2afa'JJi(p1+p3) 


Crtka 


«  A  (S)  sin(S  0i ) 

e  ■  ■■ ■■■fSTCT -  cos  (se) 

S=N,2N...  S  Bn(S) 


m=l,3 

A3n(S)  sin  (S  6i ) 


(8) 


CO 

JEL 


S=N,2N...  S  B  (S) 
'  n 


sin  (S0) 


where 


Aln(s)  =  ',S4Ak®  +  S®[3X3k4(^)  +  X]  -  (l+3k®)(Xsk®  +  vX3)j 
A3n(S)  =  iS  j-S4k®X  (^)  -S3  |V(2k®(^)  +  ^  k4)  -Ak2(MH)] 
+  £  -k3  £?)}+  X3k®(l+(£p)  v) 

+  *  (v-  7T  (l+k2 ))  ]J 

A3n(s)  =  f~  Ss  +  S4  [-2  +  ~  yv  +  3k3  (1-v))] 

+  S3  O4  ((2-v)  +  :k®  (3(l-v)  -  (^)S)j 

-  Xs  ^4-2  v  +  k®  (3~3v)j  +  1] 

+  [Xs  (l+3k®)(-~)  +  X®  ( l+k3 )  ( 3+l/k® )  O^sp)  ]  * 

J 

Bn(S)  =  js8(l+k®)  +  S6  [-2 (l+k®)  +  A®(lRk®(£|”)  +  k4  (2^))  ] 

+  S4  [X4  (6  +  3ks(2-v)  -  v®k4) 

-  X®  (2(4-v)  +  k®  (7-5v)  +  k4  3(l-w))  +  (l+k®)] 

+  S®  [X6  (4+  k®  +  2.C|~V.).  k4) 

-  3X4  (2  +  (2-v+v®)  k®) 

+  x®  (2  (2-v)  +  ■’fen),  k®  +  .3i^)k4.j  ] 

+  (l+3k®j  X4  [1+  kg-  -  2vX®  +  (1-k®)  X4]|  . 


Equations  (8)  and  (4)  thus  provide  representations  for  the  displacements 
u,  v,  w,  and  so  a  solution  to  the  problem  in  the  form  of  a  double  Fourier 
series  in  x  and  0.  It  is  readily  established  that  these  series  for 
displacements  converge  and,  in  fact,  sufficiently  rapidly  that  they  may 
be  differentiated  to  provide  convergent  series  for  the  stress  resultants 
and  couples.  To  complete  the  stress  analysis  of  the  cylinder  under 
tangential  axial  loadings  over  a  series  of  N  pads,  it  is  necessary  to 
superimpose  upon  the  stresses  due  to  the  above  self-equilibrating  load 
system  those  due  to  a  rotationally  symmetric  band  of  tangential  axial 
loads  (see  Figure  2).  We  must,  therefore,  add*  to  N  the  terms: 

ft  =  0  O^x^x-6 

X  0 

-u  (x-x  +6) 

=  ~226-2 -  xo-6  2  x  *  xq  +  6  (all  0)  (10) 

=  -N  x  +6  *  x  £  L 

o  c 

■where  N  =  26ap2,  Finally,  we  note  that  the  intensity  of  p  (x,0)  over 
the  periodic  pad  areas  then  has  a  magnitude  of  Pi+ps. 

LONGITUDINAL  MOMENT  SOLUTION 

To  account  for  the  fact  that  the  axial  load  is  not  applied  at  the 
middle  surface  of  the  shell,  we  assume  that  there  is  a  net  longitudinal 
moment  exerted  over  the  area  beneath  each  pad,  and  that  this  moment  may 
be  described  in  terms  of  a  radial  pressure  distribution  whose  intensity 
is  proportional  to  its  algebraic  distance  from  xQ: 

Pr(x,9)  =  P(6)  (x-x0)  x0-6  *  x  *  x0+5  ' 

where  p(9)  has  a  unit  value,  say  p  ,  beneath  a  pad  and  vanishes  elsewhere. 

The  solution  for  this  type  of  loading  may  be  obtained  most  expediently 
by  utilizing  Nash  and  Bridgland's[5]  basic  results  for  a  radial  line  load 
acting  over  a  portion  of  the  circumference  at  x=§.  Their  solution  for 
this  problem  may  be  written  in  the  form 


8 


w  =  2  J’in(0)  sin  (Xx) 

n=l 


v  =  E  Fsn(e)  sin  (Xx) 
n=l 


u  =  2  F3n(6)  cos  (Xx), 

n=l 


where 


r  ' 

F  (6)  =  sin  (X§)  Jp*(0)E  *(0)+  2  E  p*(S)E  *(S)cos  S0 

mn'  Eh  kaL  f  n “  3=1  11111 

m~l,3 

(  )  (12) 
Fan(e )  =  sin  (*§)  2i  E  p*(S)  Egn#(S)  sin  Bel, 

at  ksL  (  b=1  J 

p*(S)  is  the  finite  Fourier  transform  of  p(e ) , 


R  (s) 

E  *(S)  =  ^S_l 

run  '  •  " 


B  (S) 
m 


Aln(S)  =  S4(l+k8)  +  X8S8  [2+2  (l-v)k®  +  3Ftt)  k4]  +  X*[i+3k2] 

Aan(S)  =  »iS  | S8  [(l+k8)  +  k3X3  (2+  k3)]  +  Xs  [2+u+2kaX8]  (lljA) 


A3n(S)  =  A  ksS*  +  Sa  [1+3  (^)  k8  X2]  -  (l+3k8)  [u+k3Xa]  Xs 


I 


Bn(S)  =  S8  (l+k2)  +  S6  J  X2  [4 


7”3v 


k2+%l 


k4]  -2(l+ks)] 


+  S4 


+  Xs  Ss 


X4  [6+3(2-v)k3-i;2k4]  -  \s[2(4-u)+(7-5v)k2+3(l-v)k4]+(l+ka)j 
p[4+(~^)  k®+9  (~)  k4]  -  3  A2  [2+(2-w+u2)  k2]  (llffl) 
+  [2  (2-v)  +  7  £r)  +3  (~2)  k4]j 

+  (l+3k2)  X4  [1+  -  2vX®  +  (l-ks)  X4]. 


We  may  treat  this  solution  as  a  kernel  and  integrate  suitably  to  obtain 
the  desired  solution  for  the  radial  load  with  axial  variation.  For 
example,  if  the  displacement  w  of  Inference  5  is  denoted  momentarily  by 
w6(x,0,§),  then  the  displacement  w(x,0)  due  to  a  band  of  radial  pressure 
varying  linearly  about  Xq  and  within  a  band  of  length  2a6  is  given  by 


where 


v  (x,0) 


*)*  ^ 

J  (S-*o)  %  (x,8,§)  d§ 

V6 


<a 

s 

n=l 


hn  (») 


rVs 


J' 


v8 


(§-xo)sin(X§)d§f 


sin  (Xx)  (15) 


Fin(0)  sin  (X§)  =  Fln(0). 


Since 


Xq+6 

J  (S-x  )  sin  (X§)  d§  = 

x  -6 
o 


cos  (Ax  ) 
X2  0 


sin(X6)-X6  cos(^6)l , 


it  follows  lat  the  desired  solution  for  a  distribution  of  longitudinal 
moment  a  :  u  over  a  series  of  N  pads  may  be  expressed  by  Equations  11 
through  14,  provided  F  (0),  (m  =  1,2,3)  are  modified  by  being  multiplied 
by  the  quantity 


r 


f 


\ 
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-|§-  cot  (Xxo)  |sin  (X6)  -  \6  coa  (Xb)J  , 
and  provided  that  in  F^e),  p*(S)  is  given  by 

P*(S)  =  F  |p(e)j=  ^|po  sin  (se1)  JS.  ^(rn-i) (01+03 )isj 
=  ^  P0  si“  (301 )  (uj  S  =  k.U 


(16) 


(k  1,2,...) 


=  0 


S  ^  k.N  . 


i 


NUMERICAL  EXAMPLE  AND  »i3CUSSI0N 

By  way  of  illustration  of  the  precec_ng  results,  computations  have 
been  made  for  a  cylindrical  shell  with  the  parameters 

a  *  10",  L  *  35%  h  =  0.04", 

Xq  a  0.8,  E  a  107  psi  V  a  0-3 

and  having  three  areas  over  which  the  axial  load  is  exerted  with  param¬ 

eters  6  a  0.10  and  6^  a  0.025.  (Hiese  correspond  to  three  brackets,  two 
inches  in  length  and  one-half  inch  in  width,  located  eight  inches  from 
the  top  of  the  cylinder).  Figures  4a,  b,  c,and  d  present  stress  and 
moment  resultants  for  the  case  of  a  self -equilibrating  axial  load,  while 
Figures  5a,  b,  c,  and  d  present  similar  results  for  the  case  of  a  longi¬ 
tudinal  moment  M  acting  over  the  same  areas.  If  such  a  configuration  were 
used,  for  example ,  to  transmit  into  the  cylindrical  shell  a  total  load 
(force)  of  1500  pounds,  and  if  the  distance  from  the  middle  surface  of 
the  shell  to  the  application  of  the  load  were,  say,  one-quarter  inch, 
then  pj  =  500  psi  and  M  3  125  in- lb.  Uie  following  distribution  of 
stresses  on  the  inner  and  outer  surfaces  of  the  cylinder  would  result 
along  0oQ: 
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Table  I.  SURFACE  STRESSES  (PSl)  VERSUS  AXIAL  DISTANCE  (IN.) 


ax 

ax 

Aout 

ahn 

ff®out 

9 

-21,662 

+20,980 

-18,523 

+35,909 

9-1/2 

+522 

-734 

-8,491 

+12,853 

10 

-300 

-1,500 

-5,796 

+4,946 

10-1/2 

-1,482 

-1,080 

-4,025 

+2,875 

12 

-1,837 

-1,233 

-1,174 

+1,314 

In  Reference  2,  Bijlaard  presents  extensive  numerical  results  for 
localized  loading  of  cylindrical  shells  by  radial  loadB  and  by  longitudi¬ 
nal  moments,  ibis  provides  an  opportunity  to  compare  the  results  of 
Flugge's  equations  for  distributed  loads  vith  those  employed  by  Bijlaard 
for  a  wide  range  of  shell  parameters.  Furthermore,  the  comparison  can  be 
made  on  the  basis  of  actual  stresses  in  the  shell  rather  than  on  the  basis 
of  displacements.  Tables  II  and  III  typify  the  results  of  such  compari¬ 
sons.  (in  these  tables,  we  employ  Bijlaard's  notation  in  which  a  =  ^/a, 

7  =  a/h,  Pj  =:  0i,  p2  3  and  for  these  computations  Pi  =  p2 )  •  For  radial 
loads,  we  observe  that  maximum  differences  occur  when  both  p  and  7  are 
small,  and  that  as  7  increases,  Nx  and  Ng  tend  to  agree;  but  if  p  <  <1, 
results  for  Mjj.  and  Mg  by  the  two  methods  may  still  disagree  when  y  is 
quite  large. 


Table  II.  NONDIMENS IONAL  STRESS  RESULTANTS  (RADIAL  LOAD  P) 


P 

r 

aNx/P 

aN@/P 

Mx/P 

m0/p 

F. 

B. 

F. 

B. 

F. 

B. 

F. 

B. 

1/128 

5 

-1.000 

-0.9 

-0.666 

-1.1 

0.346 

0.390 

0.356 

0.430 

1/128 

50 

-9.100 

-9.0 

-9.900 

-10.0 

0.225 

0.275 

0.254 

0.330 

1/128 

300 

-55.200 

-55.0 

-57.200 

-57.0 

0.159 

0.175 

0.195 

0.250 

1/8 

5 

-0.850 

-0.8 

-0.712 

-0.9 

0.114 

0.112 

0.153 

0.155 

1/8 

50 

-6.420 

-6.5 

-3.460 

-3.5 

0.031 

0.031 

0.171 

0.070 

1/8 

300 

-25.000 

-25.0 

-5.200 

-5.3 

0.011 

0.011 

0.032 

0.031 

F.  Flugge's  equations. 

B.  Bijlaard's  equations. 
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Table  III.  NOND.TMENS IONAL  STRESS  RESULTANTS  (LONGITUDINAL  MOMENT  M) 


p 

r 

a3PNx/M 

asPNe/M 

aPMx/M 

aPMg/M 

F. 

B. 

F. 

B. 

F. 

B. 

F. 

B. 

1/32 

5 

-0.048 

-0.088 

+0.053 

-0.033 

+0.1070 

+0.108 

+0.0660 

+0.066 

1/32 

50 

-0.464 

-0.420 

-1.720 

-1.800 

+0.1010 

+0.101 

+O.O63O 

+O.O63 

1/32 

300 

-7.550 

-7.600 

-27.500 

-27.500 

+0.0780 

+0.081 

+0.0500 

+0.053 

1/8 

5 

-0.096 

-0.060 

-O.I69 

-0.250 

+0.1000 

+0.100 

+0.0600 

+0.060 

1/8 

50 

-1.860 

-1.800 

-5.630 

-5.600 

+0.0540 

+0.054 

+O.O367 

+0.036 

1/8 

300 

-10.300 

-10.300 

-18.500 

-18.500 

+O.OO98 

+0.010 

+0.0086 

+0.009 

F.  Fliigge's  equations. 

B.  Bijlaard's  equations. 
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